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Abstract — Traditional distributed source coding rarely con- 
siders the possible link between separate encoders. However, 
the broadcast nature of wireless communication in sensor 
networks provides a free gossip mechanism which can be used 
to simplify encoding/decoding and reduce transmission power. 
Using this broadcast advantage, we present a new two-encoder 
scheme which imitates the ping-pong game and has a successive 
approximation structure. For the quadratic Gaussian case, we 
prove that this scheme is successively reflnable on the {sum- 
rate, distortion pair} surface, which is characterized by the rate- 
distortion region of the distributed two-encoder source coding. 
A potential energy saving over conventional distributed coding 
is also illustrated. This ping-pong distributed coding idea can 
be extended to the multiple encoder case and provides the 
theoretical foundation for a new class of distributed image 
coding method in wireless scenarios. 

I. Introduction 

The availability of low-cost image sensor chips such as 
CMOS cameras is shifting the paradigm of sensor network 
communication from regular and small data transmissions 
to occasional and large amounts of data communications. To 
deliver images under the severe energy constraints of wireless 
sensor networks (e.g., nodes operating on batteries and a 
solar panel), multiterminal source coding has an important 
role as it can cut the rate to the theoretical lower limit. One 
possible topology for a multi-camera network is to avoid 
multi-hop routing, which can be well modeled as a typical 
distributed source coding setup where N separate cameras 
transmit correlated images to a common base station (BS). 

Looking at the typical multi-camera network setup as in 
Fig. la the cameras within the transmission range of the 



Therefore, for simplicity, we assume in this paper that all 
the overhearing nodes within the transmission range can get 
the message error-free, as long as it is received at the BS. 



emitting camera overhear the transmitted messages. We can 
take advantage of this broadcast nature, as it is given for free 
and can potentially help to simplify the coding procedure and 
reduce the transmission power. It will be shown later that 
this broadcast advantage can be used in a successive manner 
between separate encoders, which also avoids interference 
(only one node transmits at the same time). Furthermore, 
since the scale of the local camera network is usually small 
compared with their distance to the BS, it is practical to 
assume the channel capacity between neighboring cameras 
is larger than the capacity between the camera and the BS. 
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MICS, http://www.mics.orgl, a center supported by the Swiss National 
Science Foundation. 
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Fig. 1. The typical multi-camera network: (a) When Node-1 transmits 
a message to the base station, Node-2 can overhear it. (b) Setup of the 
two-encoder distributed source coding with broadcast advantage. 

To study the rate-distortion region of distributed source 
coding with broadcast advantage, we start with the simplest 
configuration of two encoders and one BS that are all within 



the transmission range of each other (see Fig. lb i: the source 



X is encoded at rate Rx without access to source Y. Based 
on the coded version of X, source Y is encoded at the 
rate Ry, which can be meanwhile regarded as a successive 
description for X. An extreme case for this scheme is when 
Rx > H(X), then it reduces to the rate-distortion problem 
of Y with side information X fully available at encoder and 
decoder. 

Successive approximation coding is another important ca- 
pability desirable for a camera network. It has been primarily 
employed in image coding, where one gets better image 
quality step by step using successive descriptions. This is par- 
ticularly useful in an energy-limited communication setup, 
since we can decide at the receiver whether a high resolution 
image is really needed after the low resolution version 
is displayed. Inspired by the ping-pong game, we extend 



Fig. lb to a new distributed source coding scheme, which has 



a successive approximation structure. As Fig. [2] shows, the 
broadcast messages act like a ping-pong ball, which is flipped 
back and fourth between the two encoders. We call such a 
scheme two-encoder Distributed Successive Approximation 
Coding using Broadcast Advantage (DiSAC2). 
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Fig. 2. Successive approximation structure of the two-encoder distributed 
source coding using broadcast advantage. 

In this paper, we present this new scheme and show that 
for the quadratic Gaussian case, it is successively refinable 
on the {sum-rate, distortion pair} surface, which is charac- 
terized by the rate-distortion region of the distributed two- 
encoder source coding (DSC2) [1]. Furthermore, we show 
that it has the possibility for rate allocation and the potential 
for energy saving. The extension of the DiSAC2 to three and 
more encoders is also briefly addressed. 

The remainder of this paper is organized as follows: 
In Sec. [II] we review the related works about multitermi- 
nal source coding and successive approximation coding. In 



Sec. Ill we present the precise setup and detailed analysis 
for a particular case of the DiS AC2 ( three-stage, with Gaus- 



sian sources and quadratic distortion). Sec. IV provides an 



extension of the new scheme to three and more encoders. 

II. Related work 
A. Multiterminal source coding 

The general multiterminal source coding problem [2] has 
been open for thirty years. The rate region for the distributed 
lossless source coding problem has been solved by Slepian 
and Wolf [3]. However, the general lossy case is not fully 
determined yet. Wyner and Ziv [4] solved a special case 
when one of two sources is entirely known at the decoder. 
Recently, Wagner et al. [1] gave the rate-distortion region 
for the two-encoder quadratic Gaussian case. Practical dis- 
tributed coding schemes [5] [6] [7] have been developed in 
recent years, however mostly restricted to sources with strong 
correlation such as 1-D measurement data, video frames, etc.. 

Distributed source coding using broadcast advantage is 
a particular form of the source coding problem with par- 



tially separated encoders. The rate-distortion region for two- 
encoder case was first addressed in [8], where coding theo- 
rems are determined for two cases: (i) one source is repro- 
duced perfectly at the receiver; (ii) one source is perfectly 
revealed to the other source. [9] and [10] further develop this 
idea to the case when the encoder in a lossless Slepian- Wolf 
setup can observe the coded data from the other encoder. 
It is proved that the admissible rate region is not enlarged. 



Fig. lb is actually the lossy case for this setup, however the 
answer for its rate-distortion region is still unknown today. 

B. Successive approximation coding 

The optimality of the successive approximation scheme for 
a single source has been studied by Equitz and Cover [11]: 
To encode a source X, a coarse description X\ with R-D 
pair (i£i, Di) is refined to a finer description X 2 with R-D 
pair (i?2,£*2)- This scheme is called successively refinable 
when any rate pair (Ri,R%) operates on the rate-distortion 
function R{D): R x = R(D{) and R x + R 2 = R(D 2 ). It is 
shown that this optimality is achieved if and only if we can 
write Xi —> X 2 ^lasa Markov chain. A Gaussian source 
with quadratic distortion is one example that is successively 
refinable on the R(D) curve. 

There is related work on successive coding for multiple 
sources. [12] proposes the sequential coding of correlated 
sources for video applications, in which the first source 
is encoded solely while the subsequent source is encoded 
based on both sources. This scheme is a weak version 
of centralized coding as it has access to both sources. 
However, it does not fully exploit the joint information due 
to the first step encoding, therefore the minimum sum-rate 
is sometimes worse than DSC2. Recently, [13] proposed 
a successive decoding scheme for the distributed source 
coding problem (no link between encoders). It is proved 
that successive decoding following a linear fusion could 
achieve the rate-distortion region of DSC2 for the quadratic 
Gaussian case. However, the final step of fusion actually 
breaks the successive decoding structure: all results have to 
be reconstructed after everything is received. 

III. Three-stage DiSAC2 

A. Setup and notations 

For the DiSAC2 scheme introduced in Fig. [2] we can 
specify the number of encoding/decoding stages. Fig. [3a] 
shows the encoding part of a three-stage setup: two sepa- 
rate encoders ENCx and ENCy cooperate using broadcast 
advantage to convey the correlated sources (X. Y) to the 
BS. The sketch for the coding procedure is as follows: 

1) At the first stage, ENCx encodes X to the codeword 
C\ without any knowledge of Y, and has a rate of 
R\. (Xi,Yi) is reconstructed at the BS after C\ is 
received. 

2) At the second stage, ENCy overhears the coded data 
C\ which is being transmitted to the BS, and it only 
transmits the refinement which fully exploits the joint 
information between the source Y and the coded data 
C\. C2 is the corresponding codeword sent in the 
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Fig. 3. Successive approximation coding schemes — encoding part: (a) Three-stage DiSAC2 for two sources, (b) Three-stage successive refinement for 
a single source. 



second stage, which has a rate of R 2 . (Xj,!^) is 
reconstructed at the BS based on (Ci,C2). 
3) Similarly, at the third stage, ENCx encodes X to the 
codeword C3 based on the coded data C\ and C2, 
which has a rate of R3. (Xa,Yi) is reconstructed at 
the BS based on (C 1; C 2 , C 3 ). 
If X = y in Fig. [3a] then it reduces to a three-stage 
successive approximation coding of a single source (see 
Fig. |3b). As we know from [11], it is successively refinable 
on the {rate, distortion} curve, as long as the Markov 
condition is satisfied. Similar results can be investigated for 
the DiSAC2 scheme. To give an insight, we specifically 



H= (0,0), 



for \p\ < 1. 



discuss the three-stage DiSAC2 (as depicted in Fig. 3a 1 with 
Gaussian sources and quadratic distortion: 

Let Af(/JL, E) be the notation for a multivariate Gaussian 
distribution with mean vector /j, and covariance matrix S. 
(X, Y) ~ Af(/i, E), where 

' 1 P 
P 1 

The corresponding distortion pairs for the three stages 
are (D Xl , D Yt ), (Dx 2 , Dy 2 ), (D x . 3 , Dy 3 ), respectively. The 
distortions are defined as Ed(X, X) where d(-,-) is the 
quadratic error measure, and similarly for Y. 

B. Coding procedure 

1) First stage: X ~ Af(0, 1) is encoded to C\ using the 
random codebook argument, with R-D pair 



R X = \\o % ~, D Xl <\. 
X\ can be decoded as: X\ = C\, thus 



(1) 



Ed(X,Xi) = D Xl - 
Since p(X) = JV(0, 1), the test channel in the first stage 

is 

p(X\C 1 )=M(C 1 ,D Xl ). 
Substituting X x = C\, 

p(X\C 1 )=M(X 1 ,D Xl ). (2) 



To decode Yi, we calculate the conditional probability 

+ CJO 

p(y|Ci) = / p(yx|d) dx 



-00 
+00 



= j p(Y\dX) ■ p(X\Ci) dX 

— OO 
+ OO 

= J p(Y\X) ■ p(X\d) dX, 



(3) 



where the third equality follows from the fact that C% is 
encoded and decoded from X (a definite function of X), 
thus p(Y\C\X) = p(Y\X). From the joint distribution of 

(X,Y), 

p(Y\X)=M{ P X.l-p 2 ). 



Substituting Q and ^ into ^ leads to 
where 

/ii = pXi, 

and 



2 7-1 2 2 

<n = d Xi p - p 



1. 



(4) 

(5) 
(6) 
(7) 



Yi can be decoded as: Y% = p,\, thus 

EdfrYt) = a\. 
To sum up, the distortion pair in the first stage is 

(D Xl ,D Yx ) = (D Xl ,al). 



(8) 



2) Second stage: C\ is known due to the broadcast 
advantage, thus according to {5}, Y-m ~Af(0, 
is encoded to C2 using the random codebook argument, with 
R-D pair 



R2 



1, < 

- log rr- 5 - 

2 & D Y ; 



D Y2 < af. 



(9) 



Y 2 can be decoded as: Y 2 — C 2 + fix, thus 
Ed(Y, Y 2 ) = Ed(Y -n^Yi-m) 
= Ed(Y-m,C 2 ) 
= D Y2 . 

Similar calculations as Sec. 



III-B.l 



where 



give 

XAl- p 2 )+D Xl Y 2P 
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D XlP 2-p2 + l ' 



(10) 
(11) 

(12) 



{D Xl( ?-p> + iy 

X 2 can be decoded as: X 2 = p 2 , thus 

Ed(X,X 2 ) = of. 
The distortion pair in the second stage is 

(D X2 ,D Y2 ) = (af,Dy a ). (13) 

3) Third stage: (Ci, C 2 ) are known due to the broadcast 
advantage, thus according to ( 10 1, X — p 2 ~ jV(0, <7§). X — 
p 2 is encoded to C 3 using the random codebook argument, 
with R-D pair 



Ra 



1, ^2 

2 



< oi 



(14) 



-X3 can be decoded as: X 3 = C3 + /i2, thus 
Ed{X, X 3 ) = Ed(X -fMt,X 3 - p 2 ) 
= Ed(X-p 2 ,C 3 ) 
= D Xa . 

Similar calculations as Sec. |III-BTT| give 

p(Y\C l C 2 C 3 )=N{pL 3 ,o 2 3 ), 

where p 3 and <? 3 are given in the Appendix. 
Y 3 can be decoded as: Y 3 — p 3 , thus 

Ed(Y,Y 3 ) = a 2 3 . 

The distortion pair in the third stage is 

(Dx 3 ,D Y3 ) = (Dx 3 ,*l). (15) 

4) Comments on coding: First, the whole coding proce- 
dure indicates that p, D Xl , Dy 2 , Dx 3 must be known both 
to encoder and decoder before coding, since (/ij, of) for i = 
1, 2, 3 are required as pre-known parameters. 

Second, special attention must be paid to the distortion 
constraints D Xl < l,Dy 2 < af,D X3 < erf, which is the 
assumptions we used for the calculation of the conditional 
probabilities. Due to the recursive structure of calculations, 
if any intermediate stage does not meet the distortion con- 
straint, negative rate would occur (which obviously is not 
allowed) which breaks the recursive chain and leads to an 
incorrect result. 

Calculations are omitted due to space, see [14] for details. 



C. Successive refinement of DiSAC2 

As we know, a Gaussian source with quadratic distortion is 
successively refinable on its {rate, distortion] curve. In this 
section, we show that a three-stage DiSAC2 with quadratic 
Gaussian case is also successively refinable on the {sum- 
rate, distortion pair} surface which is characterized by the 
rate-distortion region of the DSC2 [1] (referred to as Wagner 
Surface in this section). 

Theorem 1: For a three-stage DiSAC2 with Gaus- 
sian sources and quadratic distortion, {Ri, (D Xl , Dy^}, 
{R 1 +R 2 ,(D X2 ,Dy 2 )} and {fl x + R 2 + R 3 , (D x „ Dy 3 )} 
all achieve the {sum-rate, distortion pair} surface which is 
characterized by the rate-distortion region of the DSC2, for 
any rate triplet (R\, R 2 ,R 3 ). 

Proof: From [1], the minimum sum-rate for DSC2 is 



-Rdsc2 (D X , D\ 



- log 
2 B 
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4D x DyP 2 , 1 



1 



2D X Dy 



(16) 



At the first stage, the {sum-rate, distortion pair} is known 
from ([TJ and ([8]): 



Ri 



1, 1 

- log 

2 & D Xl 



(17) 



(D Xl) D Yl ) = (D Xl , a l). 
Combining ( [T7| ) and ( |T6] > leads to 

R% = Rdsci(D Xi , Dy^. 

At the second stage, the {sum-rate, distortion pair} is 
known from ([TJ, Q and (13) : 



i?l + i?2 



- lOff i 

2 & D Xl D Y2 



(18) 



(D X2 ,Dy 2 ) = (alDy 2 



Combining ( fT~8| > and (I61 leads to 

Ri+ R 2 = Rdsci(D X2 , Dy 2 ). 

At the third stage, the {sum-rate, distortion pair} is known 
from Q, $L4$ and ([BJ: 

1 



Ri + R 2 + R 3 



log l_L^ 

2 6 D Xl Dy 2 D X3 



(19) 



(D X3 ,D Y3 ) = [D X3 ,4). 
Combining (19 1 and ( fTo*) > leads to 

Ri + R 2 + R 3 = Rdsc2(D X3 , Dy 3 ). 

■ 

Fig. |4a| gives a visual illustration of Theorem[T[ We choose 
an initial operating point {R±, {D Xl , Dy^} at the first stage 
of DiSAC2, and then send an additional rate R 2 at the second 
stage. As the dashed curves on the gray surface suggest, 
{Ri + R 2 , (Dj,, Dy 2 )} with any R 2 operates along a one- 
dimensional curve on the Wagner Surface. 10 different initial 
points chosen give 10 curves verifying the same idea: the 
DiSAC2 is successively refinable from the first stage to the 



second stage. A similar plot (Fig. 4b 1 shows that the DiSAC2 
is also successively refinable from the second stage to the 
third stage. 




Fig. 4. Successive refinement of the DiSAC2: (a) Given a fixed Rx, the DiSAC2 with any additional rate R2 operates along a one-dimensional curve 
on the Wagner Surface (represented by the gray surface). The parameters chosen are p = 0.6, and 10 different Rx. (b) Given a fixed Ri + R2, the 
DiSAC2 with any additional rate R3 operates along a one-dimensional curve on the Wagner Surface. The parameters chosen are Ri = 0.5, p = 0.6, and 
10 different R 2 . 



D. Additional theoretical properties 

1) Feasible region of the distortion pair: As we men- 



region as: 



tioned in Sec. III-B.4 the distortion constraints Dx x < 
1)-Dy 2 < CTi,Djf 3 < it| have to be satisfied. It is shown 
in this section that given a fixed p, not all distortion pairs 
are feasible. An extreme situation is when p = 1 (Y is 
just a copy of X), the feasible region of the distortion 
pair (abbreviated as feasible region in the following) for the 
DiS AC2 is restricted to D Xl = D Yi ■ 




(a) Second stage, p = 0.3. 



(b) Second stage, p = 0.7. 




(c) Second stage, p = 0.9. (d) Third stage, p = 0.7. 

Fig. 5. Feasible region of the distortion pair. 

At the second stage of DiSAC2, we define the feasible 



{{D X21 D Y2 ) : D Xl < 1 n Dy 2 <<r?} 



(20) 



Similarly, the feasible region at the third stage of DiS AC2 
can be defined as: 

{(Dx 3 ,D Y3 ) ■. %<in D Y2 <<jjn Dx 3 < <j\) ■ 

(21) 

Fig. 5a - Fig. 5c show the feasible region at the second 



stage drawn according to (20 1. It can be seen that the feasible 
region gradually shrinks to the axis Dx 2 = Dy 2 as p 
increases, and finally ends up in the extreme situation when 
p = 1. Particularly, Dx 2 = Dy 2 <E (0,1) is always feasible, 



regardless of p. Fig. 5d illustrates the feasible region at the 
third stage drawn according to (21 1, which is shown to be 



identical with the one at the second stage (Fig. 5b 1. 

2) Rate allocation between two channels: The rate allo- 
cation between two channels is important as it is key for 
energy allocation. Recalling ([TJ, (|9j, (14i and Fig. lb we 
denote 

Rx = Ri+ R3 and R Y = R 2 

as a representation of the rate for two channels in the three- 
stage DiSAC2 scheme. 

The dashed line of Fig. [6] shows the admissible rate pair 
(Rx, Ry) of the three-stage DiSAC2, for the distortion pair 
Dx 3 = Dy 3 =0.5 and p = 0.6. It can be seen that the rate 
pair (Rx, Ry) is configurable while the sum-rate Rx +Ry 
is kept to be a constant, and there exists a particular point 
to achieve Rx = Ry- 

3) Potential energy saving: let us revisit channel capacity 
as a function of energy. For the discrete time Gaussian 
channel, the capacity [15, P.249] is 



C 



2 l0g(1 



bits per transmission, 



where E s is the energy per sample, and N s is the noise 
variance per sample. 



Fig. 6. The admissible rate pair (Rx , Ry) of the three-stage DiSAC2, 
for the distortion pair Dx a = Dy 3 =0.5 and p = 0.6. 
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Fig. 7. The energy saving (E D iSAC2 - min E DSC2 )/B DiSAC2 of the 
three-stage DiSAC2, for D x = D Y = 0.5 and different p. The DiSAC2 
is better than DSC2 in terms of transmission power. Notice that Dx\ has 
a lower bound to satisfy the distortion constraints, which varies with p. 



We consider further a band-limited channel (bandwidth 
W) with white noise (spectral density No/2), and each 
sample occupies a time interval T. In this case, the noise 
variance per sample is ^| Q -2W / 2 ^ T = No/2, Hence, the 
capacity is 



C 



1 



logfl 



2"°'' ' N a /2' 
which can be rewritten as 



bits per sample, 



,2C 



l) Nq/2 per sample. 



Due to the successive approximation structure of the 
DiSAC2, there is no interference between successive trans- 
missions. Therefore, assuming that the channels from each 
source to the base station are independent, Gaussian band- 
limited, with the same bandwidth W and noise spectral 
density No/ '2, the sum-energy for the three-stage DiSAC2 
can be modeled as 
3 



DiSAC2 



i=l 
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E 1 

8=1 
/ 1 



„2C 4 



„2fl; 
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1) No/2 



(22) 
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3 N /2, 



\D Xl D Y2 D Xs 

where the second equality follows assuming that all rates 
match the capacity for each channel usage (the Shannon 
channel coding theorem is met perfectly). 

For the DSC2 scheme, with similar assumptions as above, 
we can get 
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- 1 
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(23) 



where the second inequality follows from the fact that Rx 



Ry is constrained by the minimum sum-rate (16i. 



To compare the two schemes, we assume that the goal 
for both is to convey the source X and Y to the receiver 
while attaining a given distortion pair (Dx,Dy)- Thus, by 
plugging Dx 3 = D X ,D Y3 = D Y into (|22|, £ D iSAC2 can 
be reduced to a function of a single variable Dx 1 ■ We plot 
(£r>iSAC2 - minSDsa) /^DiSAC2 for D x = D Y = 0.5 and 
different p in Fig. [7] in which the transmission power of the 
three-stage DiSAC2 is less than DSC2. 

Since we use the non-interference model for the DSC2 
scheme, which essentially consumes less than any other inter- 
ference model (e.g., Gaussian multiple access channels [15, 
P.405]), Fig. [7] suggests that the three-stage DiSAC2 is 
potentially better than any conventional DSC2 scheme, in 
terms of transmission power. Further study of this energy 
saving effect is still in progress. 

IV. Extension to the case of many encoders 

We have addressed the DiSAC2 scheme for the two- 
encoder case in the previous sections. By similar reasoning, 
we can extend this to a general scheme with any number 
of encoders (DiSACn). Fig. [8] illustrates the setup for three 
encoders — DiSAC3: three encoders take turns to send 
messages, the BS successively decodes the messages to get 
better reconstruction of (X, Y, Z) stage by stage. The central 
principle to design the coding procedure for such a scheme 
is: each encoder should fully exploit the joint information 
between the local source, and the coded messages that have 
already been sent (including the overheard messages and 
self-generated messages). 

The minimum sum-rate for distributed source coding with 
three and more encoders is determined in [1] for a symmetric 
case: the source components are Gaussian, exchangeable, 
positively correlated, and all the target quadratic distortions 
are equal. It turns out that the scheme in Fig. [8] with three 
stages can achieve this minimum sum-rat^] It is interesting 
to study the successive refinement of this setup for this 
symmetric case and also more general conditions, which is 
the subject of on-going work [14]. 

2 Proof is omitted due to space, see [14] for details. 



% CXi,Yi,Z,) 



(X 2 , Y 2 , Zj) 



r(x 3 ,Y3,z 3 ) 



Base Station 



■ time- 



Fig. 8. Three-encoder distributed successive approximation coding using 
broadcast advantage. The straight lines represent the transmission from the 
encoder to the base station, and the wavy lines represent the same message 
that is overheard by the other two encoders. 



V. Conclusions 

We introduced a new distributed source coding scheme for 
two encoders (DiSAC2), which takes advantage of the free 
broadcast nature present in wireless networks. The coding 
procedure imitates the ping-pong game, and has a successive 
approximation structure. For the quadratic Gaussian case, 
it is proved that the three-stage DiSAC2 is successively 
refinable on the {sum-rate, distortion pair} surface, which 
is characterized by the rate-distortion region of DSC2. The 
transmission power analysis also shows a possible energy 
saving over a conventional distributed coding setup. 

We are working on the successive refinement conjecture of 
the n-stage DiSAC2, which is a generalization of Theorem [T] 
to the case n. The extension of DiSAC2 to three and more 
encoders is illustrated, and studying its properties is left for 
future work. Last but not least, the broadcast advantage is 
expected to simplify the encoding/decoding in image coding, 
therefore indicating that such a scheme has potential in 
multi-camera networks. A practical distributed image coding 
method will be the best verification of its theoretical value. 
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The explicit expression of ^ 3 and of 
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(Y 2 (p 2 - 1) (p 2 D Xl - p 2 + 1) - 

P D Y2 (X 3 (p 2 D Xl -p 2 + l)+X 1 (p 2 -I)))/ 

(p 2 D Xl (-D Y2 +p 2 -l)-(p 2 -lf). 



o-l = D Y2 (p 2 D Xl - p 2 + 1) ( (1 - p 2 ) (p 2 D X3 D Y2 + (p 2 - I)' 

p 2 D Xl {D Y2 {p 2 D X3 -p 2 + l) + {p 2 - if) J / 

[p 2 D Xl (D Y2 -p 2 + l) + (p 2 - I) 2 )' . 



